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Infectious  disease  models  play  a  key  role  in  public  health  planning. 
These  models  rely  on  accurate  estimates  of  key  transmission 
parameters  such  as  the  force  of  infection  (Fol),  which  is  the  per- 
capita  risk  of  a  susceptible  person  being  infected.  The  Fol  captures 
the  fundamental  dynamics  of  transmission  and  is  crucial  for  gauging 
control  efforts,  such  as  identifying  vaccination  targets.  Dengue  virus 
(DENV)  is  a  mosquito-borne,  multiserotype  pathogen  that  currently 
infects  million  people  a  year.  Existing  estimates  of  the  DENV 
Fol  are  inaccurate  because  they  rely  on  the  unrealistic  assumption 
that  risk  is  constant  over  time.  Dengue  models  are  thus  unreliable 
for  designing  vaccine  deployment  strategies.  Here,  we  present  to 
our  knowledge  the  first  time-varying  (daily),  serotype-specific  esti¬ 
mates  of  DENV  Fols  using  a  spline-based  fitting  procedure  designed 
to  examine  a  12-y,  longitudinal  DENV  serological  dataset  from 
Iquitos,  Peru  (11,703  individuals,  38,416  samples,  and  22,301  serotype- 
specific  DENV  infections  from  1999  to  2010).  The  yearly  DENV  Fol 
varied  markedly  across  time  and  serotypes  (0-0.33),  as  did  daily 
basic  reproductive  numbers  (0.49-4.72).  During  specific  time  periods, 
the  Fol  fluctuations  correlated  across  serotypes,  indicating  that  dif¬ 
ferent  DENV  serotypes  shared  common  transmission  drivers.  The 
marked  variation  in  transmission  intensity  that  we  detected  in¬ 
dicates  that  intervention  targets  based  on  one-time  estimates  of 
the  Fol  could  underestimate  the  level  of  effort  needed  to  prevent 
disease.  Our  description  of  dengue  virus  transmission  dynamics  is 
unprecedented  in  detail,  providing  a  basis  for  understanding  the 
persistence  of  this  rapidly  emerging  pathogen  and  improving 
disease  prevention  programs. 
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The  force  of  infection  (Fol)  describes  the  per-capita  rate  at 
which  susceptible  individuals  become  infected  with  a  patho¬ 
gen  (1,  2).  An  accurate  estimate  of  the  Fol  is  essential  for  pa¬ 
rameterizing  disease  models  (3).  It  can  be  used  to  calculate  key 
quantities  such  as  the  basic  reproductive  number  (T^o)  (2,  4)  and 
the  critical  vaccination  coverage  threshold  (p^)  of  a  pathogen  (5), 
which  are  frequently  used  to  guide  disease  control  programs  and  for 
determining  the  control  effort  required  to  eliminate  a  disease  (6). 

Dengue,  a  mosquito-borne  disease  whose  incidence  and  geo¬ 
graphic  range  have  increased  considerably  in  the  past  50  y  (7,  8), 
is  caused  by  any  of  four  related  but  antigenically  distinct  virus 
serotypes  (DENV-1,  DENV-2,  DENV-3,  and  DENV-4).  Pre¬ 
vious  estimates  of  the  Fol  for  DENV  are  few  and  uncertain 
owing  to  limitations  inherent  to  most  available  DENV  datasets, 
including  difficulty  in  specifying  when  an  individual  DENV  in¬ 
fection  occurred.  Given  the  growing  public  health  need  for  op¬ 
timal  vector  management  strategies  and  the  growing  potential 
for  deployment  of  a  dengue  vaccine  in  the  near  future  (9),  there 


is  a  pressing  need  for  accurate,  serotype-specific  estimates  of  the 
Fol  and  pc  for  DENV.  Here,  we  use  a  unique,  long-term  sero¬ 
logical  dataset  from  Iquitos,  Peru  to  provide  to  our  knowledge 
the  first  such  estimates. 

Basic  mathematical  models  of  pathogen  transmission,  such  as 
the  catalytic  model  where  the  Fol  was  initially  introduced  (1), 
make  simplifying  assumptions  about  the  parameters  governing 
transmission,  including  the  frequent  assumption  that  parameters 
do  not  vary  through  time  in  epidemiologically  important  ways 
(10).  The  assumption  that  the  Fol  is  constant  in  time  is,  however, 
inconsistent  with  current  understanding  of  DENV  epidemiology 
because  transmission  clearly  varies  seasonally  and  year  to  year 
(8,  11-14).  Resolving  the  magnitude  of  temporal  variations  in 
the  quantities  that  govern  or  summarize  transmission  requires 
(/)  adequate,  temporally  resolved  incidence  data  and  (//)  de¬ 
velopment  of  an  estimation  approach  specifically  designed  to 
use  such  a  dataset  to  compute  time-varying  quantities. 


Significance 

Using  mathematical  models  to  extend  knowledge  of  pathogen 
transmission  and  recommend  optimized  control  efforts  is  de¬ 
pendent  on  the  accuracy  of  model  parameters.  The  rate  at 
which  susceptible  individuals  become  infected  [the  force  of 
infection  (Fol)]  is  one  of  the  most  important  parameters,  but 
due  to  data  constraints  it  is  often  incorrectly  assumed  to  be 
constant  over  time.  Using  a  bespoke  method  for  a  12-y  longi¬ 
tudinal  dataset  of  serotype-specific  dengue  virus  (DENV)  in¬ 
fections,  we  estimated  time-varying,  serotype-specific  Fols  for 
all  four  DENV  serotypes.  The  Fol  varied  markedly  in  time,  which 
implies  that  DENV  transmission  dynamics  are  complex  and  are 
best  summarized  using  time-dependent  transmission  parame¬ 
ters.  Our  results  provide  more  accurate  measures  of  virus 
transmission  dynamics  and  a  basis  for  improving  selection  of 
control  and  disease  prevention  strategies. 
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Review  of  approaches  for  estimating  the  Fol  of  dengue  (3) 
shows  that,  mostly  owing  to  data  limitations,  two  approaches  have 
predominated  over  the  last  30  y  (2):  (/)  methods  intended  for  use 
with  data  collected  passively  from  clinics  and  hospitals  (clinical 
cases)  and  (ii)  methods  for  use  with  data  actively  collected  from 
age-stratified  serologic  surveys.  Data  that  are  collected  from 
verified  clinical  cases  are  temporally  resolved  but  only  capture 
clinically  apparent  infections.  Not  accounting  for  the  potentially 
large  proportion  of  undetected,  inapparent  infections  can  greatly 
diminish  the  accuracy  of  Fol  estimates  (15).  Indeed,  the  ratio  of 
DENY  infections  that  are  subclinical  (i.e.,  inapparent*)  to  those 
that  are  clinically  apparent  can  be  variable  and  often  quite  large, 
ranging  from  0.9:1-40:1  and  higher  (11,  16-18). 

Thus,  estimates  of  the  DENY  Fol  based  on  verified  clinical 
cases  (e.g.,  refs.  19-21)  are  uncertain.  In  contrast,  serological 
surveys  theoretically  capture  all  (or  most)  infections  in  a  study 
population.  The  tradeoff  is  that  the  actual  time  of  the  infection 
cannot  be  defined  from  a  single  blood  sample,  and  so  the  indi¬ 
vidual’s  age  is  used  instead.  An  important  complication  for 
multistrain  pathogens,  like  DENY,  is  that  the  infecting  strain  is 
often  not  determined  in  serological  surveys.  Investigators  in 
several  studies  estimated  the  Fol  of  DENY  using  single  sero¬ 
logical  surveys  (e.g.,  refs.  4,  13,  and  22),  but  cross-reactive  anti¬ 
bodies  obscure  identification  of  the  infecting  serotype.  Thus,  in 
most  dengue  endemic  settings  one  cannot  resolve  potentially 
important  relationships  among  serotypes  such  as  antigenic- 
dependent  enhancement  (23)  using  single  blood  specimens  from 
cross-sectional  surveys.  Owing  to  limitations  of  both  serology  and 
reported  clinical  case  data,  current  estimates  of  the  Fol  for 
DENY  are  uncertain  and  potentially  inaccurate. 

Prospective,  longitudinal  studies  generate  serial  samples  from 
the  same  individuals  over  several  years  that  can  be  used  to  val¬ 
idate  serological  results  and  derive  serotype-specific  infection 
information  (11).  A  longitudinal  study  design,  therefore,  pro¬ 
vides  data  that  are  amenable  to  estimating  the  Fol,  particularly  if 
the  Fol  changes  from  year  to  year  and  is  serotype-specific  (24). 
Because  existing  Fol  estimation  approaches  were  not  designed  to 
use  longitudinal  data,  we  developed  a  spline-based  modeling 
approach  to  analyze  a  12-y  longitudinal  serology  dataset  (11,703 
participants  and  38,416  blood  samples)  from  the  city  of  Iquitos, 
Peru  and  produced  the  first,  to  our  knowledge,  time -varying,  se¬ 
rotype-specific  Fol  estimates  for  DENY. 

Methods 

Approval  of  Experiments  Involving  Human  Subjects.  This  study  used  in¬ 
formation  from  participants  in  five  overlapping  cohorts.  Each  had  separate 
human  subjects  protocols  (see  SI  Appendix,  Table  SI  for  protocol  numbers) 
that  were  in  compliance  with  US  federal  regulations  governing  the  pro¬ 
tections  of  human  subjects.  All  protocols  received  approval  from  the  in¬ 
stitutional  review  boards  (IRBs)  of  all  participating  institutions  and  from 
a  Peruvian  Ethics  Committee  that  ensured  that  all  Peruvian  regulations 
governing  the  protection  of  human  subjects  were  followed.  Starting  in  2007, 
the  Naval  Medical  Research  Center  Detachment  (now  NAMRU-6)  formed  an 
IRB  that  is  registered  with  the  Department  of  Defense,  the  Office  Human 
Research  Protection,  and  the  Peruvian  Ethics  Committee.  In  addition,  all 
protocols  were  reviewed  and  approved  by  the  Loreto  Regional  Health  De¬ 
partment,  which  oversees  health-related  research  in  Iquitos.  In  all  instances, 
written  consent  was  provided  by  study  participants. 

Data  and  Seroconversion  Identification.  Iquitos  is  an  isolated  city  of  ^^370,000 
inhabitants  located  in  the  Amazon  basin  of  northeastern  Peru.  It  has  been 
well  described  elsewhere  (11).  All  four  DENV  serotypes  have  been  in¬ 
troduced  into  Iquitos  and  subsequently  circulated  endemically:  DENV-1 
in  1990  (25),  DENV-2  in  1995  (26),  DENV-3  in  2001  (11),  and  DENV-4  in 


*Throughout  the  manuscript,  the  terms  "inapparent"  and  "clinically  apparent"  are  used 
to  differentiate  between  asymptomatic  infections  or  infections  with  mild  symptoms  that 
do  not  result  in  detection  through  passive  case  detection  (inapparent)  and  infections 
severe  enough  that  the  individual  seeks  medical  attention  (clinically  apparent). 


2008  (27,  28).  Our  analysis  includes  12  y  (1999-2010)  of  data  from  five 
longitudinal  dengue  cohorts  involving  Iquitos  residents  >5  y  of  age  (see 
SI  Appendix,  section  SI  and  Table  SI  for  details;  Fig.  1).  In  each  cohort, 
participants  provided  blood  samples  for  testing  by  plaque  reduction 
neutralization  tests  (PRNTs)  at  6-9  mo  intervals.  In  some  circumstances, 
studies  overlapped  in  time  and  space  {SI  Appendix,  section  SI).  For  our 
analysis,  we  combined  all  cohorts  into  a  single  subsample  of  the  Iquitos 
population  with  extensive  turnover  of  individuals  throughout  the  study  pe¬ 
riod.  Some  of  the  later  cohort  studies  recruited  individuals  that  were  either 
not  yet  born  or  too  young  to  qualify  for  inclusion  in  the  first  cohort  (i.e.,  aged 
less  than  5  y  in  1999).  As  such,  to  maintain  a  comparable  subpopulation  (and 
avoid  confounding  effects  of  birth  and  a  changing  population  size  in  Iquitos), 
we  removed  individuals  from  consideration  who  were  born  after  1995.  This 
resulted  in  the  removal  of  1,465  children  from  consideration.  Thus,  our 
analyses  were  based  on  the  subpopulation  of  individuals  born  before  1995, 
which  henceforth  we  will  refer  to  as  the  sample  population. 

All  blood  samples  from  1999  through  2010  were  analyzed  for  the  presence  of 
DENV  neutralizing  antibodies  by  serotype-specific  PRNT  (11)  in  baby  hamster 
kidney  BHK21  cells  using  a  carboxymethyl  cellulose  overlay.  Samples  were 
considered  positive  when  plaques  were  reduced  70%  or  more  (PRNT70)  using 
dilutions  of  1:60,  1:80,  1:60,  and  1:40  for  DENV-1,  DENV-2,  DENV-3,  and  DENV- 
4,  respectively  (11)  {SI  Appendix,  section  S2).  Chronological  sets  of  PRNT70 
results,  coupled  with  knowledge  of  the  timing  of  DENV  serotype  introductions 
into  Iquitos,  provided  confidence  in  our  interpretation  of  serologic  results  (29). 
On  this  basis,  we  designed  an  algorithm  to  identify  infections  that  minimized 
the  probability  of  false-positive  results  by  (/)  using  serotype-specific  thresholds, 
(//■)  ignoring  all  transient  positive  results  (e.g.,  negative-postive-negative),  and 
{Hi)  eliminating  all  participants  who  had  any  instance  of  seroconversion  to 
more  than  one  serotype  in  the  same  blood  sampling  interval  (see  SI  Appendix, 
section  52  for  details).  Completely  eliminating  such  participants  from  all  anal¬ 
yses  is  quite  conservative,  but  our  method  of  identifying  seroconversions  relies 
on  investigating  the  entire  serohistory  of  an  individual.  As  such,  any  possibly 
erroneous  PRNT  casts  doubt  on  results  for  all  serotypes.  This  conservative 
approach  means  our  estimates  of  Fols  are  biased  low.  Our  approach 
allowed  us  to  identify  tertiary  and  quaternary  DENV  infections  (29)  in 
part  because  of  sequential  introductions  of  two  novel  serotypes.  Simul¬ 
taneous  virus  isolation  and  identification  from  dengue  cases  (i.e.,  case 
data)  (30)  provided  independent  validation  of  the  patterns  we  describe 
{Discussion  and  SI  Appendix,  section  52). 

Model  Description.  For  likelihood-based  inference,  we  estimated  the  pro¬ 
portion  of  the  study  population  that  had  already  been  infected  by  time  t 
[denoted  F{t)],  rather  than  estimating  the  Fol  directly.  Defining  A{t)  as  the 
Fol  at  time  t,  we  have  {SI  Appendix,  section  S3) 

t 

F(t)  =  1  -expj-  /l(L/)QfL/|.  [1] 

This  equation  implicitly  assumes  that  the  population  is  homogeneous  and 
well  mixed  (i.e.,  every  individual  is  equally  at  risk  to  be  infected  by  of  any 
infectious  individual).  This  assumption  is  imperfect  and  some  of  its  con¬ 
sequences  are  discussed  below.  The  data  were  left-censored  {SI Appendix,  Fig. 
S1A),  interval-censored  (Fig.  1),  and  right-censored  {SI  Appendix,  Fig.  SIB), 
corresponding  to  individuals  who  entered  the  study  already  seroconverted, 
seroconverted  between  two  blood  draws,  and  left  the  study  having  never 
seroconverted,  respectively.  Using  the  probability  density  function  of  the 
infection  times  of  infected  people,  f  =  dF/dt,  the  likelihood  for  the  three 
types  of  censored  data  are  given  in  SI  Appendix,  Table  S3  (where  /  denotes 
the  time  of  infection)  (31).  We  considered  each  individual's  infection  status 
to  be  independent  of  others,  which  permitted  us  to  take  the  likelihood  of 
the  data  as  the  product  of  each  individual's  likelihood  {SI  Appendix, 
section  S3). 

To  estimate  the  Fol  through  time,  we  used  a  nonparametric  spline-based 
approach  (32).  Specifically,  we  defined  a  set  of  B-splines  (33)  as  basis  functions 
for  f.  The  advantages  of  B-splines  over  other  bases,  such  as  monomials  and 
trigonometric  functions,  are  that  they  are  flexible  and  do  not  a  priori  assume 
periodicity.  Based  on  deviation  information  criterion  (DIC;  for  details  see  SI 
Appendix,  section  55  and  Table  S2),  we  identified  six  B-splines  per  year  (72 
total)  as  the  optimal  basis  for  f  for  the  four  serotypes  and  present  these 
results  in  the  main  text.  A  model  using  four  B-splines  per  year  (48  total)  also 
worked  well  {SI  Appendix,  Table  52,  section  S6,  and  Figs.  S2-S1 1).  There  were 
no  qualitative  differences  between  results  obtained  when  using  a  model  with 
four  B-splines  per  year  compared  with  a  model  with  six  B-splines  per  year. 
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Fig.  1.  Summary  of  participants  and  interval-censored  infections.  The  top  panel  shows  the  total  number  of  active  participants  across  cohort  studies  from 
1999  to  2010.  The  absence  of  a  cohort  study  from  late  2005  to  mid-2006  is  indicated  by  the  gray  shaded  region.  Remaining  panels:  After  applying  the  se¬ 
roconversion  identification  algorithm  to  the  raw  data  the  number  of  interval  censored  infections  are  plotted  against  time.  For  all,  the  midpoint  of  the  interval 
over  which  the  infection  was  censored  is  used  to  time  infections. 


We  could  not  estimate  the  function  f  before  1999  because  we  did  not 
have  data  before  that  date.  Seroconversions  that  occurred  before  1 999  were 
left-censored,  allowing  estimation  of  the  proportion  of  the  population  that 
was  exposed  before  the  beginning  of  the  study.  We  defined  k  as  the  pro¬ 
portion  of  the  population  that  had  seroconverted  before  1999: 

to 

7c  =  f{u)du,  [2] 

with  to  representing  the  beginning  of  the  study.  Using  Eq.  2,  the  likelihood 
of  both  right-censored  and  left-censored  observations  was  rewritten  using  k. 


the  chain  to  create  1,000  estimates  of  f.  This  formed  an  em  pirical  estimate  of  the 
posterior  distribution  of  f.  For  each  day,  we  then  selected  the  middle  90% 
of  the  estimates  to  form  our  Bayesian  credible  interval  (BCI)  at  that  point. 
Throughout,  in  addition  to  BCIs,  we  present  the  posterior  medians.  There 
were  two  parameters  that  a  priori  we  knew  would  have  identifiability  and 
convergence  issues:  the  parameters  corresponding  to  the  very  beginning 
and  very  end  of  the  study.  We,  therefore,  truncated  our  estimates  to  the 
region  where  our  chains  converged  {SI  Appendix,  section  S4  and  Figs.  516 
and  517).  All  analyses  were  done  with  R  (38)  and  the  R  packagefda  (39).  We 
evaluated  convergence  with  the  R  package  CODA  (40). 


Model-Fitting  Procedure.  To  estimate  the  pdf  f,  we  used  a  Markov  chain  Monte 
Carlo  (MCMC)  approach,  specifically  an  adaptive,  Metropolis-within-Gibbs  algo¬ 
rithm  (34,  35)  (for  complete  details  on  the  fitting  procedure,  see  SI  Appendix, 
section  54  and  Figs.  512-515).  We  ran  10  chains,  each  of  length  100,000,  and 
evaluated  convergence  primarily  by  monitoring  scale  reduction  factors  {SI  Ap¬ 
pendix,  Fig.  513)  (36,  37),  trace  plots  {SI  Appendix,  Fig.  514),  and  acceptance 
probability  plots  {SI  Appendix,  Fig.  515).  For  our  analysis,  we  combined  the  last 
15,000  steps  of  each  chain  and  randomly  sampled  1,000  steps  to  remove  auto¬ 
correlation.  The  parameters  were  not  independent  of  each  other,  so  to  create 
credible  regions  for  f  (and  later  the  Fol  and  1Zo)  we  used  the  1,000  sampled  steps  of 


Parameter  and  Quantity  Estimation.  With  our  estimates  of  f  we  computed  the 
proportion  of  our  study  population  (those  born  before  1995)  still  susceptible 
at  time  t,  denoted  s(t),  as 

t 

s{t)  =  ^-  +  J  f{u)duy  [3] 

to 


The  Fol,  A{t),  was  then  {SI  Appendix,  section  53) 


[4] 
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The  number  of  secondary  infections  caused  by  a  single  infectious  person  at 
anytime  t  was  the  effective  reproductive  number,  denoted  R{t).  To  estimate 
R{t)  for  the  entire  population  of  Iquitos,  we  used  our  Fol  estimates  to  cal¬ 
culate  the  fraction  of  the  entire  population  susceptible  at  any  time  t, 
denoted  Sp{t)  {SI  Appendix,  section  S3).  Using  the  estimated  mean  time 
between  successive  DENV  infections  (i.e.,  serial  interval)  of  15-17  d  (41),  we 
approximated  R{d)  on  day  d  as  the  ratio  of  the  number  of  infections  that 
occurred  between  day  d  and  day  d  +  1  and  the  average  number  of  infections 
15-18  d  in  the  future  as  follows  {SI  Appendix,  section  S3): 

rd+^8 

/  Sp{u)A{u)du 

=  - .  [5] 

3  •  j  Sp{u)A{u)du 

We  assessed  the  sensitivity  of  Uq  to  this  interval  in  SI  Appendix,  section  S5.  We 
then  used  this  approximation  of  the  effective  reproductive  number  to  calculate 
an  estimate  of  the  basic  reproductive  number,  Uq,  on  day  d  by  scaling  our  ap¬ 
proximation  by  the  fraction  of  the  entire  population  that  is  susceptible,  Sp{d): 


no{d): 


'sp{dy 


[6] 


IZq  was  used  to  calculate  the  critical  vaccination  coverage  required  to  elimi¬ 
nate  a  pathogen.  Specifically,  the  critical  vaccination  coverage  level,  pc,  sat¬ 
isfies  the  following  relationship  (5): 


Results 

Serotype-Specific  Infections.  The  final  dataset  included  38,416  blood 
samples  that  provided  serotype-specific  infection  information  for 
11,703  individuals.  Participants  provided  1  to  13  (mean  =  3.3)  se¬ 
quential  samples  an  average  of  249.5  d  apart.  We  identified  22,301 
serotype-specific  DENV  infections,  3,276  of  which  were  interval- 
censored  (Table  1).  Because  DENV-1  and  DENV-2  had  already 
circulated  in  Iquitos  for  several  years  before  the  study  began,  the 
number  of  left-censored  infections  for  those  serotypes  was  higher 
than  for  DENV-3  and  DENV-4.  Conversely,  the  number  of  in¬ 
terval-censored  infections  for  DENV-3  and  DENV-4  was  higher, 
with  DENV-3  accounting  for  more  than  half  of  all  interval-censored 
infections.  Even  though  DENV-4  was  not  detected  in  clinics  until 
2008,  there  was  still  enough  transmission  for  more  than  800  indi¬ 
viduals  that  joined  the  study  after  2008  to  seem  to  have  already  had 
neutralizing  antibody  against  DENV-4  (Table  1).  Using  the  mid 
date  between  sample  pairs  to  time  when  seroconversions  occurred, 
we  found  the  number  of  observed  infections  per  month  varied 
markedly  through  time  (Fig.  1).  Testing  conducted  in  2004  was  not 
spread  out  throughout  the  year  and  instead  occurred  at  two  times. 
As  such,  the  mid  date  of  many  individuals  occurred  within  same 
month  (July  2004)  which  does  not  necessarily  mean  that  all  of  these 
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Fig.  2.  Number  and  order  of  interval  censored  infections  by  serotype.  For 
each  serotype  the  number  of  interval-censored  infections  are  plotted  against 
year.  Note  that  for  comparison  purposes  the  scale  of  the  y  axis  is  not  the 
same  in  each  panel.  Per  individual,  these  infections  are  broken  down  by 
which  infection  they  constitute  (primary,  secondary,  tertiary,  or  quaternary). 
Because  both  DENV-1  and  DENV-2  cocirculated  before  the  beginning  of  the 
study  period,  the  majority  of  individuals  were  already  exposed  to  at  least 
one  of  these  serotypes  and  thus  most  interval-censored  infections  were  not 
primary  infections.  For  this  same  reason  (the  cocirculation  of  DENV-1  and 
DENV-2  before  1999),  there  are  considerably  fewer  DENV-1  and  DENV-2 
interval-censored  infections  than  DENV-3  and  DENV-4  interval-censored 
infections. 


individuals  were  actually  infected  in  July.  Owing  to  the  gap  between 
cohort  studies  from  late  2005  to  mid-2006  (indicated  by  the  shaded 
region  in  Fig.  1),  we  have  no  information  on  infections  that  oc¬ 
curred  during  that  period  {SI  Appendix,  section  SI).  In  total,  84.0% 
of  the  11,703  study  participants  seroconverted  to  at  least  one 
DENV  serotype  by  the  time  they  left  the  study.  Transmission 
varied  year  to  year  and  across  serotypes,  with  a  steady  increase 
in  the  number  of  postsecondary  infections  (tertiary  and  quater¬ 
nary)  later  in  the  study  (Fig.  2).  Overall,  the  majority  of  DENV-3 
and  DENV-4  infections  were  tertiary  or  quaternary  (65.2%  and 
77.4%,  respectively).  For  36.2%  of  all  individuals  who  sero¬ 
converted  to  DENV-4,  it  was  their  fourth  infection  with  a  DENV. 


Table  1.  Summary  of  censored  data  type  by  serotype 


Data  type 

DENV-1 

DENV-2 

DENV-3 

DENV-4 

Left-censored 

7,714 

7,464 

2,980 

867 

Interval-censored 

342 

408 

1,701 

825 

Right-censored 

2,527 

2,724 

5,800 

4,658 

For  each  serotype,  the  number  of  individuals  that  were  either  left-,  in¬ 
terval-,  or  right-censored.  Left-censored  individuals  entered  the  study  al¬ 
ready  seroconverted  to  the  specific  serotype.  Interval-censored  individuals 
became  infected  to  the  specific  serotype  during  their  time  in  the  study  pe¬ 
riod.  Right-censored  individuals  left  the  study  having  never  been  infected  to 
the  specific  serotype.  Note  that  none  of  the  above  columns  adds  up  to 
1 1,703  (the  total  number  of  participants).  For  some  individuals,  no  PRNT  test 
was  conclusive  for  certain  serotypes,  and  as  such  they  were  removed  from 
consideration  toward  the  calculations  concerning  that  serotype.  Addition¬ 
ally,  DENV-4  was  not  tested  for  until  2006,  and  the  smaller  number  of  cen¬ 
sored  individuals  for  DENV-4  reflects  this. 


Model  Parameter  Estimates.  The  proportion  of  the  population 
infected  before  1999,  k  {SI  Appendix,  Fig.  S18),  was  55.4%  (90% 
BCI:  53.8-57%)  for  DENV-1  {SI  Appendix,  Fig.  SIM)  and  52.7% 
(90%  BCI:  51.5-54%)  for  DENV-2  {SI  Appendix,  Fig.  SISB). 
Conversely,  k  was  essentially  0  for  DENV-3  and  DENV-4  [0.004% 
{SI  Appendix,  Fig.  S18C)  and  0.004%  {SI  Appendix,  Fig.  S18Z)), 
respectively],  which  were  introduced  later.  It  is  important  to 
note,  however,  that  DENV-4  was  not  included  in  our  PRNT 
assays  until  2006.  The  fact  that  these  estimates  were  not  exactly 
0  was  likely  an  artifact  of  the  fitting  procedure  because  no  indi¬ 
viduals  were  identified  with  a  left-  or  interval-censored  DENV-3 
or  DENV-4  infection  until  after  the  respective  introductions  of 
those  viruses.  Our  estimates  of  the  daily  probability  of  infection,/, 
showed  rough  seasonal  fluctuations  in  magnitude  across  sero¬ 
types  {SI  Appendix,  Fig.  S19).  These  estimates  were  greatest  for 
DENV-3,  particularly  in  2002-2003.  Over  the  period  of  study, 
the  susceptible  proportion  of  the  study  subpopulation,  s{t) 
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Fig.  3.  Daily  estimates  of  Fol.  For  each  serotype,  daily  estimates  of  Fol  as  well  as  the  90%  BCI  are  plotted  against  time.  The  absence  of  a  cohort  study  from 
late  2005  to  mid-2006  is  indicated  by  the  gray  shaded  region. 


{Methods),  decreased  for  all  serotypes  {SI  Appendix,  Fig.  S20). 
At  the  end  of  the  study,  s{t)  was  reflective  of  the  order  of 
serotype  introductions  into  Iquitos:  DENV-1  (first  reported  in 
1990)  =  22.9%  (90%  BCI:  21.9-23.9%),  DENV-2  (first 
reported  in  1995)  =  26.8%  (90%  BCI:  26.8-28.7%),  DENV-3 
(first  reported  in  2001)  =  32.0%  (90%  BCI:  30.9-33.2%),  and 
DENV-4  (first  reported  in  2008)  =  56.7%  (90%  BCI:  54.4- 
58.6%).  Unlike  estimates  for  the  sample  population,  the  susceptible 
proportion  of  the  entire  population  of  Iquitos,  sp{t),  was  relatively 
stable  for  DENV-1  and  DENV-2  {SI Appendix,  Fig.  S21).  However, 
susceptible  estimates  within  Iquitos  for  the  invading  serotypes 
(DEIW-3  and  DENV-4)  decreased  at  rates  similar  to  those  within 
the  sample  population  {SI  Appendix,  Fig.  S21).  The  estimated  age 
distribution  of  infections  skewed  toward  younger  individuals  the 
longer  the  serotype  circulated  within  Iquitos  {SI  Appendix,  Figs.  S22 
and  S23  and  section  S5). 

Fol.  Depending  on  year  and  serotype,  daily  Fol  estimates  ranged 
from  0  to  0.002  (Fig.  3),  with  the  hipest  estimates  being  for  DENV-3 
and  DENV-4.  Although  there  was  a  gap  between  cohorts  from 
late  2005  to  mid-2006,  we  did  identify  nonzero  point  estimates  of 
the  Fol  owing  to  the  slight  systematic  increase  in  the  proportion  of 
left-censored  individuals  that  occurred  after  that  period  compared 


with  before.  After  analyzing  the  consistency  of  estimates  across  this 
gap  {SI  Appendix,  section  S5  and  Fig.  S24),  we  found  that  the  loss  of 
data  increased  median  values  and  credible  intervals  of  estimates 
around  the  gap.  Away  from  the  gap,  estimated  Fol  values  were 
consistent  with  those  in  Fig.  3.  Further,  the  large  credible  intervals 
and  timing  of  the  estimated  peak  Fol  for  2004  (July  2004)  may  be 
an  artifact  of  the  synchronized  timing  of  blood  draws  in  2004. 

Our  longitudinal  studies  captured  the  introduction  of  a  novel 
DENV  serotype  twice.  In  late  2001/earfy  2002,  the  Fol  of  the 
recently  introduced  DENV-3  was  estimated  to  be  significantly 
nonzero,  indicating  circulation  (II).  Owing  to  the  regular  testing 
of  the  longitudinal  cohort  participants,  the  timing  of  this  increase 
was  distinguishable  from  that  of  DENV-I  and  DENV-2  (Fig.  3). 
This  pattern  was  repeated  at  the  time  of  the  introduction  of 
DENV-4  in  late  2008/earfy  2009.  In  both  instances,  the  novel 
serotype  replaced  the  existing  serotype(s).  Weekly  and  monthly 
estimates  of  the  Fol  {SI  Appendix,  Fig.  S25  A  and  B)  displayed 
similar  patterns. 

There  were  periods  when  transmission  of  multiple  serotypes 
seemed  to  synchronize.  We  computed  Spearman  rank  correlations 
on  daily  estimates  of  the  Fol  between  serotypes  and  found  that 
DENV-I,  DENV-2,  and  DENV-3  were  all  highly  correlated 
(DENV-I/DENV-2:  Pi2  =  0.77,  DENV-I/DENV-3:  Pi3=0.54, 
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Fig.  4.  Yearly  estimates  of  Fol.  For  each  serotype  yearly  estimates  of  Fol  as 
well  as  the  90%  BCI  are  plotted  against  time.  The  absence  of  a  cohort  study 
from  late  2005  to  mid-2006  does  not  preclude  the  estimation  of  yearly  Fol 
estimates  for  either  2005  or  2006,  as  evidenced  by  nonzero  Fol  estimates  for 
the  circulating  serotypes  for  both  of  those  years. 


and  DENV-2/DENV-3:  ^23=0.78).  DENV-4,  having  been  in¬ 
troduced  in  late  2008,  does  not  display  high,  if  any,  correlation  with 
the  other  serotypes  {p^^=  -  0.068,  P24  =  0.25,  and  P34  =  -  0.037). 
Fols  for  all  serotypes  were  elevated,  however,  in  late  2009/early 
2010,  with  DENV-1,  DENV-2,  and  DENV-4  reaching  their  local 
maxima  at  almost  the  same  time  in  early  2010. 

Aggregating  our  estimates  by  year  indicated  that  the  Fols  for 
DENV-3  and  DENV-4  were  highest  (Fig.  4).  For  both,  there 
were  two  years  where  the  yearly  Fol  exceeded  0.2.  The  largest 
yearly  Fol  for  any  serotype  was  the  2008  estimate  for  DENV-3 
[0.33  (90%  BCI:  0.3-0.36)].  Every  serotype  had  at  least  one 
year  with  a  yearly  Fol  that  exceeded  0.1,  corresponding  to 
seroconversion  in  10%  of  members  of  the  study  population 
that  were  still  susceptible  to  that  serotype  in  that  year.  The 
Fols  for  multiple  serotypes  were  relatively  high  in  2002,  2004, 
2008,  and  2010.  Conversely,  the  Fols  of  multiple  serotypes  were 
simultaneously  relatively  low  during  2001,  2003,  2005,  and  2006 
(relative  to  serotype-specific  values  for  the  surrounding  years). 
The  yearly  Fols  of  DENV-1,  DENV-2,  and  DENV-3  were  all 
correlated  (Spearman  rank  correlations:  DENV- 1/DEN V-2, 
P12  =  0-79;  DENV-l/DENV-3,  =  0.73;  and  DENV-2/DENV-3, 
P23  =  0-75).  Spearman  rank  correlations  with  DENV-4  were  not 
informative  because  there  were  only  three  estimated  yearly  Fols 
for  DENV-4.  As  noted  above  with  the  daily  Fol  estimates,  the 
yearly  estimates  for  all  four  serotypes  were  high  in  2010,  each 
exceeding  0.1. 

Serotype-Specific  IZq  and  Vaccination  Thresholds.  Our  estimates  of 
R{t)  and  Kq  for  each  serotype  fluctuated  temporally  (SI  Ap¬ 
pendix,  Figs.  S26  and  S27,  respectively).  Small  variations  in  the 
daily  estimates  of  resulted  in  large  variations  in  IZq  because  sp 
appeared  in  the  denominator  of  Eq.  6.  This  resulted  in  wide 
BCIs;  the  posterior  distributions  had  long  upper  tails  (SI  Ap¬ 
pendix,  Fig.  S27).  Thus,  for  the  purposes  of  comparison,  we 
plotted  estimates  of  IZq  with  the  corresponding  50%  BCI  (Fig. 
5)  to  truncate  extreme  values  on  the  upper  end  of  the  posterior 
distribution.  We  investigated  the  sensitivity  of  these  results  to 
our  definition  of  the  serial  interval  and  found  that  the  estimated 
values  were  robust  to  changes  in  this  interval  (SI  Appendix, 
section  S5  and  Figs.  S28  and  S29). 


We  did  not  estimate  a  value  of  IZq  under  1  for  DENV-1  or 
DENV-2,  except  for  a  small  portion  of  the  credible  interval  in 
a  few  instances.  The  median  values  of  IZq  for  DENV-1  ranged 
from  1.40  to  3.64.  The  median  values  for  DENV-2  were  lower 
(1.36-3.49).  Our  estimates  of  DENV-3  ranged  from  below  1  to 
4.72  by  2010.  For  DENV-4,  IZq  stayed  below  2  until  2010,  when  it 
increased  to  3.19  (Fig.  5).  Weekly  and  monthly  estimates  of  IZq 
(SI  Appendix,  Fig.  S25  C  and  D)  resulted  in  similar,  slightly  lower 
values.  There  was  a  considerable  amount  of  fine-scale  temporal 
variation  in  all  four  T^o  estimates.  Computing  the  cross-correla¬ 
tion  between  the  Fol  and  Rq  revealed  a  systematic  lag  of  '-60--80  d 
(Pearson  correlation,  SI  Appendix,  Fig.  S30v4  and  Spearman 
correlation,  SI  Appendix,  Fig.  S305).  Plotting  IZq  against  the  Fol 
(SI  Appendix,  Fig.  S31)  illustrates  that  sharp  increases  in  the  Fol 
were  preceded  by  sharp  spikes  in  IZq. 

Analogous  to  the  computation  of  yearly  Fol  estimates,  we 
computed  yearly  average  IZq  values  (SI  Appendix,  Fig.  S32)  by 
taking  the  weighted  average  of  the  daily  T^o  estimates  (weighted 
by  relative  number  of  infections).  As  with  the  daily  IZq  estimates, 
the  yearly  estimates  for  DENV-1  and  DENV-2  were  relatively 
similar,  with  DENV-1  estimates  slightly  higher.  The  yearly  IZq 
estimate  of  DENV-3  experienced  the  largest  jumps  from  one 
year  to  the  next  (2005-2006:  1.38-1.97  and  2008-2009:  1.82- 
2.61).  Averaging  over  an  entire  year  ignores  seasonality,  and  the 
highest  yearly  IZq  estimates  were  lower  than  the  highest  daily 
estimates  for  all  serotypes.  For  DENV-2  and  DENV-3,  the 
highest  yearly  IZq  estimates  occurred  in  2010  [DENV-2:  2.54 
(50%  BCI:  2.45-2.63)  and  DENV-3:  2.61  (50%  BCI:  2.55- 
2.69)],  and  the  highest  yearly  IZq  estimate  for  DENV-1  and 
DENV-4  occurred  in  2010  [DENV-1:  2.62  (50%  BCI:  2.56- 
2.68)  and  DENV-4:  1.43  (50%  BCI:  1.39-1.47)]. 

From  serotype  to  serotype  and  across  years,  values  of  IZq  and, 
in  turn,  pc,  varied.  When  DENV-3  or  DENV-4  were  first  in¬ 
troduced  (when  the  entire  study  population  was  susceptible),  the 
estimated  IZq  was  effectively  1  [DENV-3:  50%  BCI:  (0.80-1.45); 
DENV-4:  50%  BCI:  (0.76,  1.44)],  giving  pc^50%.  Using  the 
upper  bound  on  the  90%  credible  interval  for  DENV-4  at  the 
time  of  its  introduction  (2.18),  we  found pc  =  54%.  In  the  first  4  y 
of  the  study,  the  largest  median  IZq  calculated  was  3.64  for 
DENV-1  in  2001.  In  the  next  4  y,  the  largest  identified  was  3.07 
for  DENV-2  in  2004.  In  the  last  4  y,  the  largest  identified  IZq  was 
4.72  for  DENV-3  in  2010.  These  values  result  in  recommendations 
of  vaccination  coverage  of  73,  67,  and  79%  of  the  population,  re¬ 
spectively.  Using  the  largest  yearly  IZq  estimate  (2.61  for  DENV-3 
in  2009),  62%  of  the  population  would  need  to  react  to  vaccination 
with  a  protective  immune  response.  Conservatively,  following  the 
highest  estimated  TZq  overall  (4.72  for  DENV-3  in  2010),  our  results 
indicate  a  vaccine  should  be  distributed  to  79%  of  the  population. 


Discussion 

Our  results  quantify  temporal  variation  in  the  Fol  for  each  DENV 
serotype  over  a  12-y  period,  highlighting  marked  differences  in 
transmission  intensity  both  intra-  and  interannually.  We  found  that 
Fol  estimates  for  the  recently  introduced  serotypes  DENV-3  and 
DENV-4  were  higher  than  those  of  DENV-1  and  DENV-2,  which 
caused  outbreaks  before  but  not  during  the  period  investigated  > 
(1990  and  1995,  respectively).  Overall,  there  was  high  correlation  2 
between  Fols  across  serotypes.  There  were  years  of  relatively  high  S 
(e.g.,  2002,  2004,  2008,  and  2010)  and  relatively  low  (e.g.,  2001, 

2003,  2005,  and  2006)  transmission,  pointing  to  common  drivers  of 
DENV  transmission  dynamics.  Our  estimates  for  IZq  varied  from 
'-I  to  over  5,  depending  on  year  and  serotype.  During  the  years 
following  its  invasion,  DEN-4  IZq  estimates  were  less  than  those 
for  other  serotypes.  This  is  consistent  with  the  notion  that 
DENV-4  is  less  transmissible  than  the  other  three  serotypes  (42, 

43).  It  should  be  noted,  however,  that  IZq  for  DENV-4  appeared 
to  be  increasing  as  our  study  period  ended.  IZq  was  variable 
across  seasons,  warning  against  quick  estimates  for  critical 
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Fig.  5.  Daily  estimates  of  T^o-  For  each  serotype  daily  estimates  of  IZq  as  well  as  the  50%  BCI  are  plotted  against  time.  The  absence  of  a  cohort  study  from  late 
2005  to  mid-2006  is  indicated  by  the  gray  shaded  region.  The  estimates  for  both  DENV-3  and  DENV-4  are  truncated,  excluding  estimation  before  their  re¬ 
spective  introductions. 


vaccination  coverage.  Using  the  highest  estimated  Kq  values  (for 
DENV-3  from  2010),  we  conservatively  calculated pc  to  be  79%. 

Our  yearly  DENY  Fol  estimates  mostly  fell  within  the  wide 
range  of  previously  calculated  estimates  conducted  in  a  variety  of 
locations  using  an  array  of  methods  and  datasets  that  often  ag¬ 
gregated  across  DENY  serotypes.  Yearly  Fol  estimates  ranged 
from  0.07-  0.14  (19)  to  0.2-0.25  (4,  22)  in  Thailand  to  0-0.3  (13) 
in  Brazil.  Similar  to  the  Fol,  our  estimates  of  the  IZq  of  DENY 
fell  within  the  wide  range  of  previously  published  values  (1.3- 
6.3)  (44).  As  detailed  in  a  review  of  previous  Kq  estimation 
efforts  (44),  only  three  estimation  methods  used  serotype-spe¬ 
cific  data  (4,  19,  45),  and  only  one  of  those  was  based  on  sero¬ 
logical  surveys  (a  study  using  single  blood  samples  from  1,009 
children  all  collected  in  early  1980)  (4).  Our  estimates  indicate 
dengue  is  slightly  more  transmissible  than  directly  transmitted 
diseases  such  as  severe  acute  respiratory  syndrome  (7lo^2-5) 
(46)  and  influenza  (T^o  2-3)  (47)  and  less  transmissible  than  fast¬ 
spreading  diseases  such  as  measles  and  pertussis  (TIq^  12-18  and 
'-12-18,  respectively)  (48).  Although  our  estimated  critical 
vaccination  percentage  of  79%  was  high,  it  was  considerably 
lower  (and  thus  vaccination  would  be  a  more  reasonable  control 
option)  than  that  for  measles  and  pertussis  ('-92-94%). 


Our  serotype-specific  approach  revealed  synchronous  dy¬ 
namics  among  DENY  serotypes.  There  were  high  correlations  in 
both  the  daily  and  yearly  estimates  of  Fol  between  DENV-1, 
DENV-2,  and  DENV-3.  DENV-1,  DENV-2,  and  DENY-4  all 
achieved  their  local  maxima  at  essentially  the  same  time  in  2010. 
There  were  several  transmission  seasons  when,  independent  of 
the  size  of  the  serotype-specific  susceptible  pool,  there  seemed 
to  be  more  than  40  seroconversions  to  at  least  three  different 
serotypes  (Fig.  2),  specifically  in  2004,  2008,  and  2010,  even 
though  for  at  least  2004  surveillance  data  suggested  a  single  se¬ 
rotype  dominated  (30).  Even  under  stricter  schemes  for  identi¬ 
fication  of  seroconversions  (SI  Appendix,  section  S2  and  Fig. 
S33),  there  remained  periods  where  multiple  serotypes  seemed 
to  circulate  concurrently.  This  emphasizes  the  potential  for  dif¬ 
ferences  between  patterns  of  disease  (i.e.,  clinically  apparent 
infections)  and  patterns  of  infection.  A  study  identifying  the 
timing  of  serotype-specific  outbreaks  of  dengue  in  Thailand  (49) 
similarly  identified  seasonal  synchronization  across  serotypes, 
specifically  between  DENV-1,  DENY-2,  and  DENY-3.  DENY-4, 
however,  was  reported  to  be  out-of-phase.  Our  estimated  Fols 
for  DENY-1,  DENV-2,  and  DENY-4  were  at  or  close  to  their 
highest  values  at  almost  the  exact  same  time  in  2010,  indicating 
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that  the  interserotypic  immune  reactions  that  drive  patterns  of 
transmission  among  serotypes  may  vary  in  their  influence  in  dif¬ 
ferent  contexts.  We  also  noted  that  invading  serotypes  replaced 
the  existing  serotype(s)  that  had  been  previously  circulating  at 
relatively  high  levels.  Novel  serotype  invasions  are  rare  events,  in 
our  case  two  over  12  y,  which  prevented  us  from  performing 
statistical  tests  on  serotype  replacement  patterns. 

The  synchrony  between  serotypes  indicates  there  are  common 
drivers.  DENY  is  dependent  on  a  mosquito  population  to  complete 
its  transmission  cycle.  Interannual  variation  in  climate  drivers  can 
thus  impose  variation  on  transmission  dynamics  through  their  in¬ 
fluence  on  mosquito  biology  and  ecology.  As  in  many  cities  with 
endemic  dengue,  Iquitos  employs  various  mosquito  control  strate¬ 
gies  in  response  to  increases  in  dengue  cases.  We  are  currently 
investigating  relationships  between  interannual  variation  in  poten¬ 
tial  climate  drivers,  vector  control  efforts,  and  Fol  estimates. 

Although  we  computed  daily  estimates  of  TIq,  there  was  no  in¬ 
herent  reason  why  these  estimates  should  not  all  be  equal.  By  in¬ 
vestigating  which  assumptions  were  violated  to  produce  such 
temporally  fluctuating  estimates  and,  in  particular,  produce  the 
lagged  patterns  observed  between  our  estimates  of  TZq  and  the  Fol, 
it  was  possible  to  indirectly  deduce  characteristics  of  transmission 
dynamics.  For  instance,  the  repeating  pattern  of  a  sharp  spike  in  TIq 
followed  by  a  relatively  slower  increase  and  then  decrease  of  the 
Fol  was  consistent  with  violation  of  the  assumption  of  a  well-mixed, 
spatially  homogeneous  population.  We  considered  the  Iquitos  co¬ 
hort  to  be  one  subpopulation,  but  in  reality  people  (and  their  ex¬ 
posure  to  Aedes  aegypti  bites  in  locations  other  than  their  home) 
were  nonhomogeneously  distributed  across  the  city  (28,  50).  These 
heterogeneities  can  have  important  implications  for  epidemiological 
prediction  and  inference  (51).  Localized  outbreaks  occur  in  Iquitos 
(52),  and  thus  focal  outbreaks  would  best  be  scaled  by  focal  levels  of 
human  immunity.  Perhaps  if  we  incorporated  the  relatively  rapid 
depletion  of  locally  susceptible  individuals  there  would  be  a  more 
gradual  change  in  the  estimated  values  of  TIq.  Another  spatial 
heterogeneity  that  could  contribute  to  the  patterns  we  observed  was 
the  variation  in  individual  movements.  Some  people  may  have 
contributed  more  to  transmission  than  others  by  moving  about 
more,  being  bitten  by  more  mosquitoes,  being  more  infectious,  or 
some  combination  of  these  factors  (53).  Perhaps  once  an  outbreak 
was  initiated  the  pathogen  spread  to  less-transmissible  individuals, 
decreasing  the  aggregated  IZq  estimates. 

Certain  caveats  exist  with  our  approach.  First,  cross-reacting 
anti-DENV  antibodies  can  result  in  false-positive  PRNT  results 
(54,  55).  We  compensated  for  this  by  excluding  individuals  who 
seemed  to  seroconvert  to  multiple  serotypes  in  the  same  time 
period  and  by  using  an  individual’s  entire  serohistory  to  guard 
against  transient  false  positives.  This  does  not  guarantee  that  we 
completely  controlled  for  incorrectly  serotyping  an  infection 
owing  to  a  cross-reaction.  Although  there  were  periods  of  syn¬ 
chrony  between  the  serotypes,  increases,  for  example,  in  the 
number  of  DENV-3  seroconversions  were  not  systematically 
accompanied  by  an  increase  in  DENV-1  seroconversions.  As 
such,  we  concluded  that  it  is  unlikely  that  cross-reacting  anti¬ 
bodies,  rather  than  cocirculation  of  multiple  serotypes  of  DENY, 
were  responsible  for  our  results.  Second,  although  the  virus  iso¬ 
lation  and  PCR  data  confirm  certain  patterns  (e.g.,  the  timing  of 
the  DENY-3  and  DENY-4  invasions),  virus  detection  was  not  al¬ 
ways  concordant  with  the  longitudinal  serological  data.  Yirus  de¬ 
tection  in  cell  culture  and  PCR  indicated  that  one  or  at  most  two 
predominant  serotypes  at  a  time  produced  clinically  apparent 
infections.  The  patterns  of  inapparent  infections  do  not  by  defi¬ 
nition  exactly  match  those  of  apparent  infections  (15).  Considering 
the  complex  interplay  between  order  of  infection  and  severity  of 
disease,  especially  for  tertiary  and  quaternary  infections  (29,  56),  it 
was  not  surprising  that  virus  isolation/PCR  data  and  the  longitu¬ 
dinal  serologic  data  did  not  always  perfectly  agree.  Third,  we  mod¬ 
eled  each  serotype  independently,  ignoring  the  potential  effects  of 


temporary  cross-protection.  Following  a  DENY  infection,  an  in¬ 
dividual  has  temporary  immunity  to  heterologous  DENY  infection 
(57),  which  may  have  affected  our  estimates  (44).  Correcting  for 
cross-protection  would  result  in  a  systematic  increase  in  our  Fol 
estimates  because  temporarily  immune  individuals  would  be  re¬ 
moved  from  the  heterologous  susceptible  pool  in  the  denominator  of 
Eq.  4.  Likewise,  incorporating  death  would  increase  our  estimates. 

Although  they  require  considerable  effort,  time,  and  resour¬ 
ces,  longitudinal  studies  provide  valuable  detailed  information 
on  pathogen  transmission,  especially  when  rates  of  asymptomatic 
infection  are  high,  like  they  are  for  DENY  (16).  Because  most 
estimation  attempts  are  not  based  on  longitudinal  data,  the 
methods  to  use  such  detail  are  not  well  developed.  As  noted 
earlier,  the  use  of  likelihood-based  fitting  of  smooth  functions  for 
the  Fol  has  been  developed  for  data  from  a  single  serological 
survey  of  individuals  within  the  study  population  (2).  By  com¬ 
bining  a  spline-based  approach  with  equations  previously  derived 
in  the  field  of  reliability  (31)  we  developed  an  estimation  method 
that  incorporates  the  detailed  information  provided  by  a  longitu¬ 
dinal  study  and  allows  for  temporal  flexibility  in  Fol  estimates. 
Our  method  was  designed  to  work  with  a  particular  dataset  but 
could  be  adapted  to  other  longitudinal  studies  for  dengue  and 
other  infectious  diseases. 

Our  analysis  and  the  interpretation  of  serological  data  were  fa¬ 
cilitated  by  the  two  novel  DENY  introductions  that  took  place 
during  the  period  of  study.  In  both  cases,  the  novel  serotype  initially 
displaced  the  preexisting  serotypes.  After  this  initial  phase,  however, 
neither  of  the  novel  serotypes  seemed  to  interfere  with  the  others.  If 
Iquitos  begins  to  sustain  simultaneous  transmission  of  all  four 
serotypes,  the  maintenance  of  longitudinal  cohort  studies  will  pro¬ 
vide  valuable  data  for  confirming  the  patterns  we  identified  and/or 
reveal  further  complexity  in  the  interplay  among  serotypes. 


Conclusions 

Beyond  clear  intra-annual,  seasonal  variation,  the  observed 
temporal  variation  in  epidemiological  parameters  (especially 
from  year  to  year)  implies  that  the  transmission  dynamics  of 
DENY  in  Iquitos  are  complex  and  cannot  be  summarized  with 
synoptic  data  or  assumptions  of  time-independent  transmission 
parameters.  In  addition  to  informing  dengue  prevention  strat¬ 
egies,  connecting  the  variation  of  these  estimates  to  other  pro¬ 
cesses,  such  as  measures  of  entomological  risk  (i.e.,  mosquito 
abundance)  and  climatic  variation,  will  inform  local  control 
strategies.  The  minimal  vaccine  coverage  required  to  effectively 
control  dengue  varies  by  serotype.  Given  a  potentially  limited 
number  of  doses,  efficiently  distributing  a  vaccine  within  and 
between  communities  can  only  be  optimized  when  these  varia¬ 
tions  are  understood  and  taken  into  account.  Because  vaccines 
may  not  have  perfect  efficacy  (9),  the  incorporation  of  accurate, 
serotype-specific  estimates  of  critical  transmission  parameters 
will  be  crucial  for  selecting  delivery  strategies  and  determining 
the  optimal  mix  of  vaccination  and  vector  control  for  sustainable 
prevention  of  dengue. 
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